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Abstract: In this note we have tried to determine how the existence of a local 
entropy current with non-negative divergence constrains the second order transport 
coefficients of an uncharged fluid, following the procedure described in [1]. Just 
on symmetry ground the stress tensor of an uncharged fluid can have 15 transport 
coefficients at second order in derivative expansion. The condition of entropy-increase 
gives five relations among these 15 coefficients. So finally the relativistic stress tensor 
of an uncharged fluid can have 10 independent transport coefficients at second order. 
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1 Introduction 

Fluid dynamics is an effective description of near equilibrium physics. It captures 
the dynamics of locally equilibriated systems in which the parameters of equilibrium 
vary slowly compared to relaxation length scale. When, for instance, microscopic 
dynamics is well described by kinetic theory, the Boltzman equation reduces to the 
equations of fluid dynamics on length scales that are large compared to the molecular 
mean free path. 
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The variables of fluid mechanics are the local values of the parameters that 
characterize fluid equilibrium; in the simplest context these are the fluid temperature, 
chemical potentials and velocity. The equations of fluid dynamics are simply the 
conservation of stress-tensor and all other charged currents, once the stress tensor 
and currents are expressed in terms of equilibriation parameters. The formulas that 
express the stress tensor and charge currents as functions of fluid variables are known 
as constitutive relations. As fluid dynamics is a long wavelength effective description, 
it is meaningful only to present constitutive relations in an expansion in derivatives 
of the fluid variables. 

For any given fluid, a microscopic computation of constitutive relations starting 
from a microscopic description of the system is often an impossibly difficult task. 
In the usual spirit of effective field theory, this task is, moreover, extraneous to the 
study of fluid dynamics. An autonomous 'theory' of fluid dynamics addresses the 
following question: what is the most general form of the constitutive relations that 
could possibly arise in the fluid description of any consistent system (to any given 
order in the derivative expansion). 

The requirement of symmetry restricts the form of the constitutive relations. For 
example, the stress tensor is a tensor. At any given order in the derivative expansion, 
there exists only a finite number of (onshell inequivalent) tensor structures one can 
build out of fluid variables and their derivatives. The most general expression for the 
stress tensor is clearly given by a linear combination of these inequivalent tensors, 
where the coefficients in this expansion are arbitrary functions of the scalar fluid 
variables (temperatures and chemical potentials in the simplest situations). While 
the requirements of symmetry are certainly necessary, they are not sufficient. There 
is atleast one additional constraint on allowed constitutive relations: that they are 
consistent with a local form of the second law of thermodynamics [3]. In other 
words any given constitutive relation must be accompanied by an entropy current, 
also constructed out of fluid variables. The entropy current must have the property 
that its divergence is positive for every conceivable fluid flow allowed by the fluid 
equations with the specified constitutive relations 1 . 

It is a quite remarkable fact that the requirement of the existence of a posi- 
tive divergence entropy current constrains the allowed constitutive relations of fluid 
dynamics in a quite dramatic manner; as we will see in some detail below, this re- 
quirement reduces the number of free parameters (or more precisely free functions) 
allowed in constitutive relations. 

In this note we will work the restrictions imposed by this requirment on the 
constitutive relations of an uncharged relativistic fluid in 3 + 1 spacetime dimensions 
at second order in the derivative expansion. In this system the variables of fluid 

1 However, the existence of a local entropy current with positive divergence is really a heuristic, 
with as yet no solid basis in thermodynamics or QFT 
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dynamics are simply the temperature T and the fluid four velocity u M . As we will see 
in some detail below, symmetry considerations allow a 15 parameter worth of consti- 
tutive relations for the stress tensor, where every parameter is an arbitrary function 
of the temperature. It was already noted by Romatschke [1] that the requirement 
of entropy increase imposes atleast two relations between these 15 functions. In this 
note we generalize Romatschke's analysis and demonstrate that a complete study of 
the requirements of positivity of the entropy current imposes 5 relations on the 15 
coefficients described above. The set of all 2nd order constitutive relations consis- 
tent with the positivity of entropy increase is parameterized by ten functions of the 
temperature. But unlike the first order case we did not find any inequalities among 
the second order transport coefficients. 

We now give a detailed presentation of our final results, i.e. a precise character- 
ization of the 10 parameter set of allowed constitutive relations at second order in 
the derivative expansion for uncharged relativistic fluids. Let the fluid stress tensory 
be given by 



■ in/ 



We work in the so called Landau frame which imposes the transversality condition 



w M IV = 



In this frame the most general allowed form for upto second order in the deriva- 
tive expansion is given by 2 
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2 Our convention for Riemann tensor is the following 
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where 



u^ 1 = The normalised four velocity of the fluid 
pnv _ giiv _|_ u n u v _ p ro j ec ^ or perpendicular to 

6 = V.w = Expansion, a M = (u.V)i^ 
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For any tensor A 
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It turns out that 'entropy-positivity' does not impose any constraint on r, Ao, Ai, A 2 , d, £i 
and ^2- The rest of the eight second order transport coefficients satisfy the following 
5 relations. 
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So finally there are 10 independent transport coefficients at second order for some 
uncharged fluid. 

As we have explained above, unless the equations (1.3) are satisfied, the fluid 
dynamics equations do not have a positive divergence entropy current. When the 
equations (1.3) are satisfied the fluid equations are compatible with the existence 
of a positive divergence entropy current, but this current is not unique. Consider a 
fluid with a particular constitutive relation that obeys the equations (1.3). It turns 
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out that any such fluid has a 7 parameter (7 arbitrary functions of temperature) set 
of positive divergence entropy currents. 



J^ 1 1 upto 2nd order 

= su» + V u [4i(u"VT - u u WT)} + V u (A 2 W) 

+ A 3 (rt - l -g^R^j u u + + ^ [6VT - P Qb (V^)(V a T)] (L4) 

+ {B^ 2 + B 2 2 + B 3 a 2 )u^ + 5 4 [(Vs) V + 2sOV M s] 

Our results above apply to an arbitrary theory, i.e. a theory with an arbitrary 
equation of state. The specialization of our results to the special case of conformal 
fluids turns out to be trivial, as we now explain. As was explained in ([4, 5]), the 
requirement of Weyl invariance forces 10 linear combinations of the 15 symmetry 
allowed coefficients to vanish; specifically 

K\ = 2k, 2 = k, t = 3Ao, A4 = 

(1.5) 

Ci = C2 = 6 = 6 = 6 = 

Moreover, the the temperature dependence of the remaining five coefficients (those 
that are allowed to take arbitrary values consistent with Weyl invariance) is deter- 
mined by dimensional analysis; All of them are just linearly proportional to temper- 
ature. 

It turns out that this linear dependence on temperarture, the conformal equation 
of state and (1.5) reduce all of the equations (1.3) to trivial identities of the form 
= 0. In other words, the requirement of positivity invariance of the entropy current 
does not impose any equations on the five transport coefficients (allowed by conformal 
symmetry) . 

The coefficients for a conformally covariant entropy current are given by the 
following expressions. 

A\(T) = a 1} A 2 (T) = a 2 , A 3 (T) = ^-T 

B^T) = bxT, B 2 (T) = ^T, B 3 (T) = b 3 T (1.6) 



where all Oi and b{ are constants. 

Therefore the conformally covariant entropy current has four independent coefficients 
(di, a 2 , bi and 63) when expanded upto second order in derivatives. 

Let us end this introduction with a description of our motivations in undertaking 
the computations described in this note. Our first motivation is practical. Theoret- 
ical reconstructions of the RHIC and LHC heavy ion experiments often model the 
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expansion of the hot dense deconfined plasma by the equations of fluid dynamics 
including second order corrections. Given that we have not been able, from first 
principles, to compute the fluid description of QCD it seems of interest to parame- 
terize the most general set of allowed equations, as we have done (to 2nd order) in 
this note. 

However our main motivation in undertaking the computations described in this 
note are structural. At zero temperature, the equations of motion of physical systems 
are strongly constrained by the requirement that they follow from the extremization 
of an action. On the other hand the tradiational formulation of the equations of 
fluid dynamics is at the level of the equations of motion. It seems likely that the 
equations of fluid dynamics inherit constraints that are the analogue of the zero 
temperature requirement of being obtained from an action. One possible source of 
such constraints, as described in this note, stems from the requirment that our system 
admit an entropy current of positive divergence. There are other potential sources of 
constraints, for instance the requirement that correlation functions computed from 
fluid dynamcs have certain symmetry properties that can be derived, on general 
grounds, in quantum field theories (see [6]). It appears to us to be of interest to find 
a complete 'theory' of fluid dynamics; a formalism that ennumerates all consistency 
conditions on the equations of fluid dynamics. Such a formalism would take the 
place of the zero temperature requirement that the equations of motion follow from 
an action. The computations presented in this note may be thought of as a small 
first step towards this larger goal. 

2 Brief Summary of our Procedure 

In the rest of this note we proceed to determine the most general second order con- 
stitutive relations and second order entropy current consistent with positivity of the 
divergence of the current. In order to do this we first list out the most general sym- 
metry allowed entropy current upto third order in the derivative expansion (onshell 
equivalent currents are not treated as distinct). We then perform a brute force com- 
putation of the divergence of this entropy current, keeping all relevant terms (see 
below for an explaination of which terms are relevant) to fourth order in the deriva- 
tive expansion. We then use the equations of motion (including constitutive terms 
upto second order in derivatives) to rewrite our final answer entirely as a function 
of derivatives of velocity and temperature and the background metric that are inde- 
pendent of each other (i.e. are not related to each other by the equations of motion). 
We then work out the conditions on the entropy current and constitutive relations 
that ensure the positivity of this divergence for arbitrary values of the independent 
fluid derivatives, finally obtaining (1.3) and (1.4) 

In order to implement the programme outlined in the paragraph above, in this 
note we proceed in the following order. In section 3 below we classify and ennumerate 
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the onshell independent derivatives of fluid fields (upto fourth order in derivatives). 
We also enumerate t all the products of these derivative fields with net derivative 
number < 4, organizing our ennumeration in representations of the local SO (3) that 
leaves the fluid velocity fixed. 

In section 4 below we then proceed to ennumerate the most general entropy 
current upto third order in the derivative expansion. We then compute the divergence 
of this entropy current and determine several constraints on the entropy current that 
follow from the requirement that its divergence is positive definite. 

In section 5 below we then ennumerate the constraints on constitutive relations, 
upto second order, that follow from the requirement of positivity of divergence of the 
entropy current. 

We end this brief section by listing our conventions. Throughout this note we 
work in the Landau gauge. In this gauge the velocity at any point is defined as 
the unique time-like eigenvector of the stress tensor, normalised so that u^u^ = — 1. 
In other words, by definition 

T;u u = -eu, (2.1) 

The quantity e is taken by definition to be energy density of our fluid. All other fluid 
thermodynamical quantities (like the temperature or pressure) are obtained from e 
using thermodynamics and equation of state. Equation of state expresses the energy 
density e as a function of some thermodynamic parameter like entropy density and 
it can vary from system to system . In this note we shall keep it arbitrary. Once 
e(s) is known, the temperature (T) and the pressure (P) can be determined in the 
following way. 

Both of the above relations directly follow from equilibrium thermodynamics. 



3 Classification of fluid data 

3 In this section we present a partial listing of the onshell independent derivatives 
of fluid fields (T and -u M ), at any given point x, upto fourth order in the derivative 
expansion. We organize these derivatives (which we will often refer to below as inde- 
pendent data) according their transformation properties under the SO (3) rotational 
group that leaves u^{x) invariant. 

In order to explain what we mean let us consider a listing of independent data 
at first order in the derivative expansion. Before accounting for onshell equivalences 
we have 16 independent pieces of first derivative data; (the four derivatives of tem- 
perature and the four derivatives of each of the three independent velocities). These 
16 pieces of data transform, under the local SO(3), as two scalars, two vectors, a 

3 This section has been worked out in collaboration with Shiraz Minwalla and Tarun Sharma. 
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pseudo vector and a traceless symmetric tensor (i.e. the 5) of 50(3) (see the second 
column of Table 1 for details). However these 16 pieces of data are not all inde- 
pendent. The four perfect fluid equations of motion may be used to solve for four 
of these fluid derivatives in terms of the other 12. As the four fluid equations can 
be decomposed into a vector and a scalar of 50(3) (see the third column of Table 
1) it follows that the independent data consists of one vector, one scalar, a pseudo 
vector and a traceless symmetric tensor (see the fourth column of Table 1). The 
choice of the independent scalar and vector piece of data is arbitrary; we could take 
our independent data to be either of the vectors and either of the scalars listed in 
the second column of Table 1. In the fourth column of Table 1 we have made one 
particular choice of the independent data that we will employ in much of our note. 
Occasionally we will find it more convenient to use another choice of independent 
data; we will explicitly point this out when this is the case. 

In this note we will require that the production of entropy is positive for an 
arbitrary fluid flow on an arbitrary curved manifold. As explained in [7] this require- 
ment yields constraints for the constitutive relations of fluids even in the flat space. 
In order to implement the constraint described above, we will find it necessary to 
list the data assocaited with local background metric curvatures in addition to the 
data associated with fluid flows. All curvature invariants formed from a background 
metric are given in terms of (contractions of) the Riemann tensor and its deriva- 
tives. It is important to recall that, in addition to certain symmetry properties, the 
Reimann tensor also obeys a Bianchi type identity. The independent derivatives of 
the Reimann tensor should be counted modulo the Bianchi identity and its deriva- 
tives. In analogy with the counting problem for fluid data listed above, we will regard 
the set of all derivatives of the Reimann tensor (with all symmetries imposed) as raw 
data, and the Bianchi identities and its derivatives as 'equations of motion'. We will 
then list the indpendent pieces of data in curvature and derivatives by subtracting 
equations of motion from raw data, just as described in the previous paragraph. 

3.1 Independent Data 

With no further ado we simply proceed to list the (relevant parts of) fluid and 
curvature data at various orders in the derivative expansion. 

At first order in derivatives we have only fluid data. They are listed below in 



Here for any tensor A^, the symbol Ai^ means the symmetric traceless part 
of it, projected in the direction perpendicular to u^. 



For example, if we expand this notation, the shear tensor a^ v has the following 
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Table 1. Data at 1st order in derivative 





Before imposing eom 


Eoms 


Independent data 


Scalars (1) 


(u.V)T, (V.u) 


u v V^ v = 


= (V.u) 


Vectors (1) 




P£V v T va = 


= (w.VK 


Pseudo- vectors (1) 


u u e u ^ Xu V x u a 




^ = t^e^V^ 


Tensors (1) 






0> = V^iv) 



definition 



pa po 



V a «f> + VfcWa 







2 """3 

At second order we have curvature data along with fluid data. The curvature 
data is given by the 20 independent components of the Reimannn curvature subject 
to the identities 



R^iua/3 


-R 


Rfj,vaf3 


Ra/3 




-- 



-R 



The 20 independent components may be decomposed into 50(3) representations as 
in Table 2 



Table 2. I2 type curvature data 



Scalars (2) 


R — R^ 1 fiui 

R 00 = u^vTR^ 

= U^U R fiau 


Vectors(l) 


P^ a R ab u b 


Tensors (2) 


R^v) 

where = u^R^p 


Pseudo-tensor 


UbR( l J' Cd £v)cdqU q 



Second order fluid data is tabulated in Table 3 

Third order fluid data that transforms in the scalar, vector and pseudo vector 
representations is tabulated in Table 4 (we will never need 3rd order data in the 5, 
7 and 9 representations, and so do not bother to tabulate these below). 
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Table 3. I2 type fluid data 






Before imposing eom 


Eoms 


Independent data 


Scalars (1) 


(n, X7\Gi V72T" 


u u {u.v jv ^1 p — u, 
V M V,T^ = 




Vectors (2) 


P^(u.V)a U) P^V 2 u y , 

p" v v v e, p^v u {u.v)t 


P£(u.V)V„T ua = 0, 

UeP&VbVvT™ = 


P' ia V a 0, 
P»V b a ab 


Pseudo-vectors (0) 


(u.V)P 


Uvt^V & a TZ = 




Tensors (1) 


P» a P» b {u.V)(T ab , V <M V,)T 


V (M V a T" = 


P» a P vb (u.V)a ab 


Pseudo-tensors (1) 


^ {jj.lv) 






Spin-3 (1) 









Table 4. 73 type fluid data 





Before imposing eom 


Eoms 


Independent data 


Scalars (1) 


(u.v) 2 e, v 2 e, 

(w.V) 3 T, (w.V)V 2 T 


u v {u.V) 2 V ^ = 0, 
(w.V)V M V,T^ = 0, 
^V 2 V^T^ = 


(u.v) 2 e 


Vectors (1) 


P^ a (ii.V) 3 M a 
P^ a (n.V)V a e 

P^ a V 2 V a T 
P^ a (u.V) 2 V a T 
P^ a (u.V)V 2 u a 


P^W a V M T^ = 
P£V 2 V u T ua = 

^VaV^T^ = 

P^{u.VfV v T va = 


P^(u.V)V a G 


Pseudo- vectors (1) 


P^ a (w.V) 2 / a 
P^ a V 2 / a 


Ui/C /«'^( u .V)V a V a 2^ = 


P^ a V 2 / a 



The third order curvature data consists of derivatives of the Reimann curvature 
constrained by Bianchi identity 

e abcd \/ b p t/aj3cd = 

In Table 5 we list the independent curvature data that transforms in the scalar, 
vector and pseudo vector representations (again we will not need and so do not list 
the remaining representations) 

Finally, fourth order scalar data (all we will need at fourth order), both fluid as 
well as curvature, is tabuated in Table 6. 
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Table 5. I3 type curvature data 





Before imposing eom 


Eoms 


Independent data 


Scalars (2) 


(u.V)R 
(u.V)Rqo 


u^ aa pe abcd V b R af} cd = 


(u.V)Rqo 
(u.V)R 


Vectors (3) 


P^ a R 00 

P» a V a R 
P^ v F va 
P%V v R va 
p^u\u.V)R ab 


UaU^^VbR^cd = 

UaU^e^e^VbR^cd = 


P^ a R 00 

P^ v R va 
P^u h {u.V)R ab 


Pseudo- vectors (1) 


uPu a e abcd V b R% cd 
v?u a e» abc V b R pc 


u a u a e abcd V b R a ^ cd = 


u P Ua€ »abc VbRpc 



Table 6. I4 type scalars 





Before imposing eom 


Eoms 


Independent data 


Fluid data (1) 




(w.V) 3 KV m T^) = 


(u.V) 3 




(w.V)V 2 6 


(w.V)V 2 (u v V^ u ) = 






(w.V) 4 T 


(m.V) 2 V m V v T^ = 






(u.V) 2 V 2 T 


V 2 V M V„T^ = 






V 2 (V 2 T) 






Curvature data (4) 


(u.V) 2 R Q0 
(uV) 2 R 
V 2 R 

V 2 i? 00 
u a {u.V)V b R ab 
V a V b R ab 
V Q V fe F ab 


u a u u e^V^ abcd V b R^ cd = 

v"w v "^v 6 ^ d = 

^^•V^V^^ = 


(w.V) 2 ^ 
(u.V) 2 Roo 
V 2 R 

V 2 i?oo 



3.2 Composite Expressions 

In the sequel we will sometimes need to list for example, all 3rd order vectors. In 
addition to expressions constructed out of the independent 3rd order data, listed in 
the previous subsection, the set of all 3rd order vectors includes expressions cubic 
in first order data, and expressions formed out of the product of one first order 
and one second order piece of data. We will refer to expressions constituted out 
of products of independent data as composite expressions. Composite expressions 
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formed out of independent data are easily ennumerated and decomposed into 50(3) 
representations using Clebsh Gordan decompositions (taking care to acccount for 
symmetry properties when we multiply two or three copies of the same data). 

In order to ease the process of reference to composite expressions in the rest of 
the note we now adopt the following terminology. Independent data at m th order in 
the derivative expansion is referred to data of the type I m . A composite expression 
that consists of a product of three first order pieces of data is referred to as an 
expression of the type C\^\. Composite expressions that consist of the product of a 
first order and 3rd order piece of data are called expressions of the form 0^3. The 
generalization of our notation to other forms of composite data is obvious. 

In the sequel we will need to list only those composite expressions that transform 
in the vector (in order to list the most general entropy current) or the scalar (in order 
to list the most general terms in its divergence). In the rest of this subsection we 
present a listing of those vector and scalar composite expressions that will be needed 
below. 



Table 7. C\ \ type expressions 



Scalars (4) 


6 2 , a 2 , u 2 , a 2 


Vectors (3) 




Tensors (5) 


©CT/iy, vtpPav)-, ^l^av), ^^au), &<M C M 



Here = P^P vb [ v ° tt '- v '"° ] 



Table 8. C± t 2 type expressions independent of the curvature 



Scalars (4) 


e( w .v)e, (v m t)v 2 < (v^TXti.vxwr), ^vvt 


Vectors (11) 


P;eV 2 M", P*Q(u.V)V»T, P*(V 2 T)V»T, u^V 2 ^ 
w^-VJVT, a^W, a^u.V^T, P;(V a T)(V^V a T) 
(V a uvK 6 , P;a Qb V a V 6 < P«u ab Wu ab 



Table 9. C\ 2 type expressions involving a curvature 



Scalars (5) 


F a b& ah ', R a b& ab , u a a b R ab , OR, QRoo 


Vectors (9) 


P^F^, PfavRr, 0^00, a"i2, P«u a eR a » 
u a R ab a b ^ u a R a » bc uj bc , P«u a R ab » c a bc , u a R ab u b » 
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Table 10. Cx,i ; i type expressions 



Scalars (7) 


6 3 , a 2 6, u 2 Q, a 2 Q 
a^tyr"". a^tfa^, u^a^ 


Vectors (10) 


aV, wV, 6V, aV, Qa^a u 
Quj^a u , o» a o ab a\ u» a u ab a\ u» a a ab a\ a» a u ab a b 



4 Entropy current 

In this section we will derive constraints on constitutive relations, at second order 
in the derivative expansion, from the requirement of positivity of divergence of any 
entropy current that reduces to su^ (where s is the entropy density) in equilibrium. 

We will first explain in very broad terms how we proceed. 

The entropy current takes the form 

J" = J£ + > (4.1) 

where 

In general J M has terms of all orders in the derivative expansion, but for the purposes 
of this note we will find it sufficient to truncate J^ 1 to terms of third order or lower in 
derivatives. To start with we allow J to be given by the most general possible form 
consistent with symmetries. We then compute the divergence of J M and reexpress 
the final result entirely in terms of independent data of fourth or lower order in 
derivatives. The last step (reexpressing the divergence of J in terms of independent 
data) uses the equations of fluid dynamics, and so the constitutive relations. Our 
final expression is a polynomial in the (finite number of) pieces of data of fourth 
or lower order in derivatives. We then demand that the resultant polynomial is 
positive definite (or can be made so by the addition of terms higher than fourth order 
in the derivative expansion) as a function of its arguments. This rather stringent 
requirement turns out to yield several constraints on the form of both the entropy 
current at second (and third) order as well as constitutive relations at second order 
in the derivative expansion. 

4.1 Entropy current in equilibrium and 1st order correction 

In order to set the stage for our discussion we first recall how the requirement of 
positivity of the entropy current constrains constitutive relations at first order in the 
derivative expansion [7]. We first recall that thermodynamics and the fluid dynamical 
equations may be used to demonstrate that 

i ( ©n^\ 

V^ = --(n^ + ^J (4.2) 
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It follows in particular from (4.2) that entropy is conserved in perfect fluid dynamics 
(i.e. when IP" vanishes). It also follows that the divergence of the most general 
entropy current (4.1) only contains terms of second or higher order in the derivative 
expansion. Let us now examine the constraints from the requirement of positivity 
of these second order pieces. For this purpose we need to study the most general 
entropy current at first order in derivatives. Imposing the requirement of invariance 
under partiy, the most general (onshell inequivalent) family of first order entropy 
currents is given by 

J" = ^equilibrium + ^ + W (4-3) 

(we have used here that at first order in the derivative expansion we have one piece of 
scalar data, which may be chosen as 0, and one piece of vector data, which may be 
chosen as (P). We now proceed to compute the divergence of (4.3) and use the perfect 
fluid equations to rexpress the result in terms of independent data. The resultant 
expression is the sum of a quadratic form in first order data and a linear form in 2nd 
order scalar data. As derived in [7], the final expression for this divergence is given 
as 



V^t/^lupto 2nd order 
1 

= ~ T 



( B\ Q 2 (44) 
+ Q(u. V)a + (a. V)/3 + la + | J — + f3 (a 2 + u 2 ) 



+ (a + (3)(u.V)0 + fiRoo 

Here both the first and the second line have terms quadratic in 1st order data. 
The last line contains the terms which are linear in second order data. There are 
three independent 2nd order scalars ((w.V)O, R, Roo) as given in the classification 
in section 3. Only two of these three scalars appear in (4.4). Since these two terms 
are linear in fluid variable, they can have any sign and to ensure positivity of the 
divergence their coefficients (both a and (3) must be set to zero. This implies that at 
1st order no correction can be added to the entropy current which is consistent with 
the positivity requirement. Then in the RHS of (4.4) only the first line will give a 
non-zero contribution. To evaluate the first line we need the first order corrections to 
the constitutive relation. At first order the most general correction to the constitutive 
relation (stress tensor in Landau gauge) will involve the single on-shell independent 
1st order scalar which we have chosen to be G and single on-shell independent tensor 

| upto 1st order V^^v Q^P^iv 

where r] and £ are shear and bulk viscosity respectively. 
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Therefore finally 

V M J M | upto 2nd order = 7f (V°* + C© 2 ) ( 4 - 5 ) 

Hence to have a positive definite divergence one requires that 

t] > o, c > o 

The main point to note in the above equation (4.5) is that it involves only 
two of the four first order on-shell independent data as listed in section 3. The 
squares of the independent vector a M and the pseudo-vector P do not appear in 
equation (4.5). Because of this fact any term in the divergence which is of the form 
(o M x I 2 or I 3 type vector) or (Z M x I 2 or I 3 type pseudo-vector) can never be made 
positive-definite. 

4.2 General constraints on second and the third order corrections 

In general J M can be written as 



J" 



where (3, is an arbitrary combination of ith order on-shell independent scalars and 
93f is a combination of ith order vectors. In the previous subsection we have seen 
that to constrain the first order transport coefficients rj and C we need to determine 
only the first order correction to the entropy current (i.e. only @i and 23^ and both 
of them finally turn out to be zero). But to constrain the second order transport 
coefficients we need to go till the third corrections to the entropy current. The reason 
is the following. 

Suppose the divergence of the most general entropy current has two terms of the 
form 

Va J*t = Ax 2 + Bxy = Ax 2 ( 1 + ^ 

\ Ax 

where x and y are two on-shell independent fluid data and A and B are some functions 
of temperature, which in general will depend on the coefficients appearing in the 
entropy current or transport coefficients. 

In this schematic expression of divergence since x and y are two independent 
fluid data, locally the ratio ^ can take any negative value, larger or smaller than 1 
in magnitude and the positivity constraint will depend on whether y 2 term is present 
or not in the final expression of the divergence. In the absence of a y 2 piece, the 
coefficient B has to be set to zero and the coefficient A to some non-negative number. 

But this argument does not require x and y to be of same order in derivative 
expansion. Even when x is of first order in derivative and y is of second order, 
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the ratio can be of order 1 for some particular fluid configuration where x is 
accidentally small enough to be comparable to y at a given point. In such cases to 
see whether y 2 term is present or not, we need to compute the divergence till fourth 
order. This is why we have to compute the divergence till fourth order even if we 
want to constrain just the second order transport coefficients. 

In fact the constraints on transport coefficients will involve situation where x 
and y are necessarily of different orders. For example, B will contain some second 
order transport coefficients only when x is of first order (as we will see below that 
x has to be equal to a^ v or 0) and y is of second order in derivatives. It will turn 
out that most of the equalities among the coefficients will follow from this sort of 
argument. 

Below we schematically list all the constraints we need to impose on the third 
and fourth order pieces of the divergence in order to ensure its positivity. 

• The coefficient of any term (appearing in third or fourth order piece of the 
divergence) which contains more than one factors of a^ v or or at least one 
factor of 00^ will not have any constraint from positivity as long as r\ and 
£ are non-zero and are of order one. This is because whenever such third 
or fourth order terms are non-zero, the second order piece of the divergence is 
also non-zero and positive-definite and will always dominate these terms within 
derivative expansion. These terms can never make the divergence negative. 
Therefore while calculating the the third and fourth order divergence we shall 
ignore all these terms. 

• One needs to do sixth order analysis to constrain the coefficients of any term 
(appearing in the fourth order divergence of the entropy current) which is of the 
form [piiv x Some third order tensor) or (0 x Some third order scalar). Such 
terms generically will have contributions from third order transport coefficients. 
Since we are interested only upto second order transport coefficients we shall 
ignore all such terms while calculating the fourth order divergence. 

• The coefficients of all the terms which contain a single I2, I3 or 1^ type scalar 
(at second, third and fourth order respectively) have to be set to zero. This 
is because locally all these terms are linear in fluid variables and therefore can 
have any sign. 

• In the third order and fourth order piece of the divergence, the coefficients 
of all the terms which are of the form (o M x I2 or / 3 type vector) or (Z M x 

^3 type pseudo-vector) have to be set to zero. Since there is no on-shell in- 
dependent I 2 type pseudo-vector, there could not be any term of the form 
(If! x J 2 type pseudo- vector). 
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• At this stage the terms appearing in the third order piece of the divergence will 
be of the following form. 

1. o^y x (J 2 or Ci t i type tensors) 

2. 9 x (J 2 or Ci 5 i type scalars) 

All these terms will involve the second order transport coefficients. 

The relevant terms appearing at the fourth order (where all the terms involving 
(Tpv and G are ignored) will be of the following form. 

1. A quadratic form involving independent 1% type data 

2. A quartic form involving a M and 00^ 

3. Terms linear in I2 type data and quadratic in a M and/or 

Therefore when 77 7^ and ( 7^ the relevant part of the divergence calculated 
upto fourth order is schematically given by 



Divergence = — 

+ <yy,v x (h or Ci t \ type tensors) + x (I 2 or C^i type scalars) 
+ A quadratic form involving independent I2 type data 
+ Terms linear in I 2 type data and quadratic in a M and/or u^ v 
+ A quartic form involving and 

(4.6) 

where in the second line all the C\ t i type tensors involving and all the C\ t \ 
type scalars involving are ignored. 

Now we can shift by a combination of I2 or C\ t \ type tensors such that 
the term linear in a^ v appearing in the second line gets absorbed. This shift 
will generate fourth order terms structurally similar to the terms appearing in 
third, fourth and fifth line of the above equation. One can see that all these 
newly generated terms together will necessarily be negative definite. Similar 
shift has to be done to absorb the terms linear in to the first line of (4.6). 

One can do similar shifts in I 2 type data to absorb the terms appearing in the 
fourth line of equation (4.6) into terms appearing in the third and fifth line with 
I2 data replaced by the shifted one. At this stage the schematic expression of 
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the divergence will take the following form. 

T) (shifted of + C (shifted 6) 2 
Divergence = — 

+ A quadratic form involving shifted I2 type data 
+ A quartic form involving and uj^ iu 

• The positive defmiteness of the divergence finally will imply the positivity of 
the quadratic and the quartic form appearing in the second and the third line 
of (4.7). 

Such condition will generically give some inequalities among the coefficients. 
However suppose by explicit computation one finds that for some particular negative 
definite term generated by the shift there is no term present in the third or fifth 
line of equation (4.6) to compensate. Then this will imply that the coefficient of the 
corresponding linear term (the source for generating this particular negative-definite 
term through shift) in the second line or fourth line has to be set to zero. This will 
give strict equalities among the coefficients. 

It will turn out that all of the constraints on the 2nd order transport coefficients 
will arise from this sort of argument. 

In explicit calculation we will see that in the quadratic form involving the I2 type 
data there will not be any term proportional to Rl , R 2 , F^F^ and R llv R a bP tia P vb . 
This will imply that the coefficients of all the terms linear in Rqq, R, F^ y and 
R ilv P lia P ub have to be zero. It turns out that once we set these linear terms to 
zero, the quartic form mentioned in the last line of equation (4.6) also vanishes. 

The vanishing of these terms at fourth order gives the final constraint on the 
transport coefficients. In the explicit computation we will see that there are eight 
terms (6a 2 , 6 / 2 , a^a", a llv lH v o^Bf, a^F^, RQ and i? o©) in the third order 
divergence which are linear in the set of fluid and curvature data mentioned above 
and also involve eight independent transport coefficients. So setting the coefficient 
of these linear terms to zero, we can express the eight transport coefficients in terms 
of the coefficients appearing in the second order entropy current. It will turn out 
only three of the entropy current coefficients appear in these expressions. Eliminat- 
ing these three coefficients we get the final five relations among the 15 transport 
coefficients as presented in (1.3). 

Once all these relations are imposed on the divergence, one is left with a quadratic 
form involving only I2 type data . To ensure that this quadratic form is positive- 
definite the coefficients appearing in the second and third order entropy current as 
well as the transport coefficients have to satisfy some inequalities. But in this case, 
at least upto this order the entropy current coefficients can not be eliminated from 
the relations. Therefore unlike the first order transport coefficients the second order 
ones do not satisfy any inequalities within themselves. 
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4.3 Implementing the general rules at second order 

At second order we have to determine ©2 and QJ2 such that the divergence calculated 
upto third order in derivative expansion is non-negative. Here we shall follow the 
general procedure described in the previous subsection. 

We shall express ©2 and QJ^ i n terms of the on-shell independent second order 
scalars and vectors respectively. ©2 will have 7 coefficients, three multiplying the 
three independent I 2 type scalars and rest of four multiplying the four C\ t \ type 
scalars. will also have 6 coefficients, three multiplying the three I 2 type vectors 
and the rest multiplying the three C\ \ type vectors. 

So before imposing any constraint the entropy current at second order contains 
total 13 coefficients, each of which is an arbitrary function of temperature. 

We shall write this most general 13 parameter entropy current in the following 
form . 



Here s is the entropy density and all the coefficients At and the Bi are the arbitrary 
functions of temperature. 

Now we shall argue that (4.8) is actually the most general 13 parameter entropy 
current. By equations of motion one can show that the only I 2 type vector appearing 
in the first term is P Ma V a O, the second term contains a linear combination of all the 
three independent I 2 type vectors and in the third term the only I 2 type vector that 
appears is P fia R a i ) u b . Therefore the first three terms together take care of the all 
the three I 2 type vectors. Terms multiplying A4, A 5 and A 6 are the three I 2 type 
scalars. 

By equation of motion B4 term is equal to a linear combination of 2 M At , a 2 u^ 
and Oa'\ Similarly B^ term is a particular linear combination of Q 2 ^, a 2 ^, Oo^, 
c^cr^ and a v u^ u . Therefore all the C\ \ type scalars and vectors appear in (4.8) with 
distinct coefficients. 

Next we shall compute the divergence of this 13 parameter entropy current con- 
structed in (4.8). We have to set the coefficients of all the I3 type on-shell independent 
terms to zero. Since there are total 3 independent I3 type scalars, it can impose at 
most three relations among the coefficients appearing in the second order entropy 
current. Next we have to isolate all the type terms which are of the form of a M 
times a I 2 type vector and set their coefficients to zero. Since there are total three 
second order I 2 type vectors this condition also can impose at most three constraints. 

• The divergence of the first two terms in (4.8) vanish identically. 



•/^second order =V„ [il(«"VT - U»WT)} + V v (A 2 T<J") 




(4.8) 
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• The divergence of the third term (the term with coefficient A 3 ) does not produce 
any ^3 type scalar. The divergence of this term is explicitly computed in (B.2). 

• The three independent I 3 type scalars (it a it 6 V a Vi,0, u.VR, u.Vi?oo) are pro- 
duced from the three terms multiplying coefficients A4, A5 and A§ respectively. 
Therefore to maintain positivity A±, A§ and Aq have to be set to zero. 

• The divergence of the terms multiplying Bi, B 2 , B 3 and B4 do not 

produce any term of the form a M times an I 2 type vector. The divergence of 
these terms are explicitly computed in (B.5), (B.6), and (B.7) respectively. 

• The divergence of the terms with coefficients B e and B 7 produce the two terms 
di/V^cr^ 1 " and a^V^Q respectively whose net coefficient should be zero to ensure 
the positivity of the divergence. 

Since these are the only places where these terms are produced, Bq and Bj are 
set to zero. 

• Both the terms multiplying B§ and A3 produce the third possible term of the 
form a M times an I 2 type vector which is a^R^Uy (see (B.2) and (B.3)). The 
net coefficient is (4* + ^ - 

Therefore positivity implies 

dB 5 _ A 3 dA 3 
dT T dT 

After imposing all these constraints the final form of the second order entropy 
current is given as 

J^ 1 1 second order 

=V„ [Ax(u^V v T - u u WT)} + V„ (A 2 Tu^) 

+ a 3 (V - IfTR) u u + + ^ [evr - p a& (v^)(v a T)] (49) 

+ (Bxu 2 + B 2 Q 2 + B 3 a 2 )u^ + B 4 [(Vs) V J + 2s6V M s] 

4.4 Implementing the general rules at third order 

• First we have to write ©3 and QJ3 in terms of the on-shell independent data. 

• The coefficients of all the 7 4 type data appearing at the fourth order divergence 
have to be set to zero. An I 4 type term in the fourth order divergence can occur 
only when the derivative acts on the I 3 type terms of the third order entropy 
current. Therefore this condition will constrain the coefficients of the I 3 type 
terms. 
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Now there are 3 I 3 type independent third order scalars and 4 I 3 type indepen- 
dent third order vectors and there are total 5 Z4 type fourth order scalars. This 
means one can have at least 2 distinct coefficients multiplying ^3 type terms in 
a third order entropy current with positive definite divergence. 

It turns out that after imposing this constraint there are exactly two coefficients 
left and the terms multiplying them can be chosen in such a way that their 
divergence vanish identically. So these two terms will not contribute to any 
further constraint. 

• The number of free coefficients that can multiply the type data in third 
order entropy current is quite large. There can be total 9 free coefficients in 
63 and 20 in QJ^. 

These coefficients will be constrained by the fact that any terms of the form 
a M times a l 3 type vector or times a I 3 type pseudo-vector have to be set to 
zero. 

Since there are total 4 I3 type vectors and total 2 I 3 pseudo-vectors, it can 
produce at most 6 constraints, reducing the number of free coefficients to 23. 

• Since there is no other general constraint to simplify the form of the entropy 
current at this stage we have to calculate the divergence. 

But we will not attempt to calculate the full divergence. Instead we shall 
calculate only those terms which can impact the constraints on the second 
order transport coefficients. 

For this purpose in the divergence we can ignore all those terms which are 
multiplied by or a^ y . 

Also to simplify we shall try to write the terms in the form of (V ) where 
A^ u is an anti-symmetric tensor, so that their divergence vanish identically. We 
could not do it for all the independent terms, but we try to apply this trick for 
as many terms as possible. 

Now we shall explicitly construct the required part of the third order entropy 
current piece by piece. 

The part multiplying the ^3 type terms can be written as 

J 1 3rd order /I3 type 

=V„ [Fi(/(«.V)V'T - u»(u.V)VT)\ + V„ [P 2 (u»R%u - u^R u e u e )} 
+ [P 3 {u.V) 3 T + P A (u.V)R + P 5 {u.V)R 00 \ u» 
+ P» a [P e V a R 00 + P 7 V a R] 

In the first and the third term the two independent ^3 type fluid data are chosen 
to be (w.V) 2 V /i T and (-u.V) 3 T. The independent ^3 type curvature data are chosen 
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from the list given in section 3. The second term (with the coefficient P2) contains 
the independent vector P^ a u b {u.y)R a b- 

Here all the terms in the third and the fourth line produce independent I4, type 
scalars at fourth order and therefore they are all set to zero. The divergence of the 
first two terms vanish identically. 

So finally this part of the entropy current has only two terms and both the terms 
have zero divergence. 

1 3rd ordcr// 3 type ('4 11") 

=V„ [P 1 {u v {u.V)V i T - u"{u.V)V v T)] + V v [P 2 (u u R%u d - u»R v e u 6 )] 

The part multiplying the type terms have total 29 coefficients to begin with. 
We shall try to write them in a way so that the computation becomes simpler. In 
table (11) and table (12) we have listed each of the independent C li2 type fluid 
data and then the independent combination through which this data has entered 
the entropy current. In table (15) we have listed the relevant Ci,i,i type scalars and 
vectors and also their coefficients in the entropy current. In all these cases, to begin 
with the coefficients are some unspecified functions of temperature. 



Table 11. type Scalars (Fluid data) 



Scalars as listed before 


Combination that enters entropy current 


0(u.V)G 


Qi [9(«.V)9] u" 


ct^(u.V)ct^ 


Q2 [a ab (u.V)a ab ] u» 


<VV a <7<^ 


V„ [Q 3 (u»o av - U V G^) Oa] 




V„ [Q 4 {uW - uW) oJ 



We shall start our analysis by computing the divergence of the terms appearing 
in table (11) and (12). 

• The divergence of the terms with coefficients (Qi, i = 3, • • • ,7) vanish identi- 
cally. 

• It turns out that the term with coefficient Qg is the only term which produces 
a M times a third order J 3 type vector. Therefore Q$ has to be set to zero. 

• Similarly if we analyse only the third order entropy current, the term with 
coefficient Qg is the only term that produces times a third order I3 type 
pseudo-vector in the fourth order divergence. However a similar term is also 
produced when the divergence of the £>i term in the second order entropy 
current is computed upto fourth order. 

V M [EiwV] = Biw 2 6 + [(u.^B^u 2 + 2Biu; a6 (ii.VH a (4.12) 
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Table 12. d j2 type Vectors (Fluid data) 



Vectors as listed before 


Combination that enters entropy current 


ev^e 


V„ \Q 5 ®(u»g av - uV)a-l 




V„ [Q 6 Qu^} 






a a (V^V a T) 


g 8 a a (V^V a T) 


0J ab V»uj ab 




a^(M.V)6 


Qio [a"(ii.V)G - u"(a.V)6] 




^v,(g n e) 








Ql3 ^VfltT 8 " 


P?<7 ab W {c a ab ^ 




P» c a ab (V a u bc - ?fV k u kc ) 


g 15 P^a afe (V a uj bc - ?fV k uj kc ) 



Since there is no second order on-shell pseudo vector, ui ba {u.V)uj ba must contain 
a third order pseudo- vector times Z M . Then in the final fourth order divergence 
the total coefficient of such term (ie. the term proportional to third order 
pseudo- vector times 1^) will be a linear combination of Bi and Qg, which should 
be set to zero. 

But to simplify the calculation instead we shall introduce a third order shift 
in the l B\ term' of the second order entropy current and will consider the 
following term \B\uj 2 u^ — ^-u^ b VJ^ b ] . The divergence of the shifted l B\ 
term' no longer contains the terms proportional to third order pseudo-vector 
times I p. (The releveant part for the divergence of the shifted l Bi term is 
computed in (B.5). Once this shift is done, Qg also has to be set to zero, 
since now this is the only term which produces Z M times a third order /3 type 
pseudo- vector. 

For the rest of the 8 terms we have to compute the divergence explicitly. However 
we are interested in those terms in the fourth order divergence which does not have 
any explicit factor of or a^ v . This simplifies the calculation. For example, in 
the divergence of the term with coefficient Qi, the only contribution which will be 

relevant for our purpose is (o\ [(-u.V)©] 2 ^ . 

• The relevant part of the divergence of the last two terms with coefficients Q14 
and Qi5 are the following. 
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Qib P vc [VX1 



P, 



V K u kc 
(4.13) 



Here in the first line we get a term proportional to (spin-3) 2 and in the second 
line we get a term proportional to (pseudo-tensor) 2 . It will turn out that 
such terms cannot occur in any other place. Positivity of the divergence will 
be satisfied if both Qu and Q15 are positive. Therefore these terms will not 
produce any constraint on the second order transport coefficients. 

The other five terms where the relevant parts are easy to calculate are the 
following. 



v„ [Qi u"e(u.v)e] Qi[(u.v)e] 2 

V M [Q 2 u»a ab (uV)a ab ] => Q 2 [(u.V)(J ab }[(u.V)a 

V M [Qi 2 a£V v e] Qi 2 [V, t <][V^0] 

V M [Q 13 <V,a^] Q13 [V M <] [Vetr 9 "] 



(4.14) 



The relevant part in the divergence of the terms with coefficients Q10 and Qn 
are more complicated. 

The divergence of the 'Qn-term' is given by the following expression. 



v M [u^v, (Q n e)] 

= Q n [V^H [V„6] + [V^l [V v Qu] O 



- Qn [V"0] 

- Qhw 2 (m.V)0 



(4.15) 



where in the last line we have used the identity (A. 5) and ignored the terms 
proportional to 6 and a^ v . 
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The divergence of the 'Qio-term' is given by 



VJ Qio K(^.V)e-^(a.V)6] 



dQi 



-T\^3\ a 2 ( u .V)0 + s ^ ] 0(a.V)6 



+ Qio [(V.a)(«.V)0 + a"(V/)(V a 8) - (V 6 0)(u.V)a b - 0(a.V)0] 



r- 



Qio a 2 ( M .V)6 



+ g 



10 



a; 2 ( M .V)6 + [(w.V)6] 2 + R 00 (u.V)Q - P ab (V a Q)(V b Q) 



(4.16) 

To express the divergence in the chosen basis of independent data we have used 
the identities (A.l), (A. 3) and (A.4). 

Here also in the final expression we have ignored the terms proportional to 
and a^ u . 



Table 13. C\ t 2 type Scalars (Curvature data) 



Scalars as listed before 


Combination that enters entropy current 


F ab <J ab 


Pi (F ab a ab ) u» 


Rab<? a 


P 2 (R ab a ab ) u» 


QR 


P 3 u^QR 


QRoo 


p 4 w^ePoo 


u a a b R ab 


P 5 u»(u a a b R ab ) 



Next we shall compute the divergence of the curvature type data appearing in 
table (13) and (14). 

• It turns out that the divergence of the terms multiplying P 5 , P 6 and P7 are the 
only terms which produce the terms of the form a M times an independent third 
order ^3 type curvature vector. Therefore we have to set P5, Pq and P7 to zero. 

• Similarly the divergence of the term multiplying P 8 is the only place where Z M 
times a third order I 3 type curvature pseudo-vector is produced. Therefore P 8 
should also be set to zero. 

The divergence of the remaining 10 terms have to be computed. Here also we shall 
ignore any term that are multiplied by an explicit factor of G or a^ v . 
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Table 14. C\ % % type Vectors (Curvature data) 



Vectors as listed before 


Combination that enters entropy current 


a^R 


P 6 a^R 


a^Roo 


P 7 a^R 00 


u a R ab u bfl 


P 8 u a R ab u b ^ 




P 9 (R^ - \g^R) a u 


a u F^ u 


P 10 {F^a u - R 00 a^ + u»u a a b R ab ) 


u a R a » bc u bc 


P n (u a R a ^u a(3 + 2co^u a R aa ) 


u a BR afl 


P 12 u a eR a » 


u a R ab a b » 


P 13 u a R ab o b » 


u a R ab ^a bc 


P u u a R ab ^a bc 



First we shall determine the relevant part of the divergence of the terms with 
the coefficients Pi, i = 1,2,3,4 and 12, 13, 14. These are easy to calculate 

V M [Pi {F ab a ab ) u"] => Pi F ab (u.V)a ab 
V M [P 2 {R ab a ab ) J] P 2 R ab (u.V)a ab 



v„[p 3 u»eR}^ p 3 p( M .v)e 

V M [P 4 n^ePoo] Pa Poo(n.V)0 

V,, [P12 P£u a QR ab ] P 12 u a R ab P»V,e 

V M [Pis n a P ab a^] P 13 u a R ab P b V^a 



(4.17) 



V„ [P14 u a P^P a6 ^V 6c ] P 14 u a P»R ab » c V,a bc 

Now we shall compute the divergence of the difficult terms multiplying the coef- 
ficients Pg, Pio and Pu respectively. We first analyse the situation where in a given 
basis all the fluid data are locally zero upto the required order and only the curvature 
data are turned on. For such configurations it will turn out that only these three 
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terms can produce non-zero divergence. They are given by the following expressions 4 . 

R.Roo 



P 9 (V-^R)a b 



Pa 



R ab F, 



ab 



(4.18) 



V M [P 10 (P^a, - P 00 a M + ^ a a 6 P a6 )] = P 10 [P ab P ab - R 2 Q0 ] (4.19) 



V M [P n (uaR^u^ + 2u» a u a R aa )] 

1 



Pn 



2u a u°R nr P ca R, 



A" 



-A 



A 



(4.20) 



Up RP°**pgpup\ 



where 

Now from these three equations (4.18), (4.19) and (4.20) we can conclude the 
following. 

• The last term in the RHS of (4.20) contains only the (pseudotensor) 2 and a 
(vector) 2 , but cannot produce any term proportional to Pqo or FabF ab . There- 
fore to have positivity of the divergence for all values of F 2 and Pq we must 
set Pio to zero . 

• Once Pio is set to zero, there are no terms in the final expressions of divergence 
that contain \R^ V A , (P^y)) , R? or Pq . Therefore in the full divergence the 
coefficients of all the terms linear in these four independent data must be zero. 
To satisfy this condition we have to set the following coefficients to zero. 

1. Pi = as it is the total coefficient of the term F ab (u.'V)a a b. 

2. P 2 = as it is the total coefficient of the term R ab (u.'V)o~ a i,. 

3. P3 = as it is the total coefficient of the term R(u.V) < d. 

4- Qio + P4 = as it is the total coefficient of the term P 00 (w.V)6. 
5. Pg = as it is the total coefficient of the term R.R 0Q . 



4 In this computation, apart from the explicit curvature there is one more source for producing 
the curvature terms. These are arising because we want to write the final answer for the fluid 
data in a given basis as chosen in section 3. For example while computing the left hand side of 
equation (4.18), we shall get a term like i?(V.o). However, our basis of independent second order 
fluid data contains a single scalar (u.V)0. Therefore we have to express (V.a) in terms (u.V)O 
before setting the fluid data to zero. In this process we shall generate a curvature term Rqq as 
calculated in equation (A. 3). Similar techniques have been used to compute the divergence of the 
other two terms multiplying Pio and Pn. 
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Ci,i,i type data 

In the fourth order divergence we are not interested in any terms that are 
multiplied by © or o^. Therefore in the third order entropy current we did 
not need to consider the C\^\ type terms which contains more than one factor 
of 0, a^ w or both. 

Here we are listing only those terms which we shall require for our analysis. 
Table 15. C\ \ \ type data (Only the relevant ones) 



Scalars 


Vectors 




K b u» a a ab a b 




K 6 a» a ou ab a b 


K 3 u» {a b a c a bc ) 


K 7 u» a a b ® 


K, (u ab a b u ca ) 


K 8 wV 




K 9 a 2 a" 




K 10 u^co ab a b 



Now we shall compute relevant part in the divergence of the each of the relevant 
term. It will turn out that analysing the explicit expression of the divergence 
we can further set some of the non-zero coefficients to zero. 



Here also the relevant parts are easy to calculate for the terms with the 
coefficients K iy i = 1,2, 3, 4, 5, 6 and 7 . 

V„ [Kx u^u 2 ®] K x u 2 (u.V)Q 

V M [K 2 wV6] => K 2 a 2 (u.V)& 

V M [K 3 u» (a b a c a bc )] K 3 a b a c (u.V)a bc 

V M [K 4 U " (u ab a b c u ca )] K 4 u ab [(u.V)a b ] 

V M [K 5 u» a a ab a b ] K 5 u» a [V^ ab ] a 1 

V M [K 7 u» b a b Q] => K 7 [u^a 6 V M 6] 

The divergence of the terms with coefficients K s , K 9 and K w are compli- 
cated. These are given by the following expressions. 



b 



(4-21) 
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'Ka-term' 



V, [K 8 wV] 

-- [(a.V)A 8 ] u 2 + 2A 8 cu^(a.V)u UIM + A 8 w 2 (V.a) 



dK 8 
dT 



aV + 2K 8 w""(a.V)w 



(4.22) 



+ K 8 u 2 [u 2 + (u.V)e + i? o] 



In the last step we have kept only the relevant terms and used (A.l) and 
(A. 3) for simplification. 

Kg-term 



V M [K 9 aV] 
= a 2 (a.V)K 9 + A 9 a 2 (V.a) + 2A 9 a»a v W^a u 



- T 



JKg 

It 



+ Kg a 4 + K q a 1 



(u.v)e + n 



on 



(4.23) 



+ 2K 9 a»a u [(u.V)<v + + w^w^j 



In the last step we have used equations (A.l), (A. 3) and (A. 7). 
K w -term 

V M [A 10 u;^ afe a b ] 



A" 



10 



(V^Rfca" + u^ a a 6 (V M u; a6 ) + a;^ a6 (V M a 6 



+ (V M A 10 K a u; afc a 6 



A 10 J (a M a;^ a6 a b ) + A 10 co^a b (V^ ab ) 



(4.24) 



+ iTio w"^"" [(u.V)<v + + u^a/J 



3 



In the last step we have used relevant part of equations (A.l), (A. 6) and 
(A.8). 

Now as explained before, all the terms that are linear in Rq Q , R, R^ and F^ 
should be set to zero. This will imply the following for the G\^\ part of the 
entropy current. 

1. A'g = as it is the total coefficient of the term u 2 Rqq. 

2. Kg = as it is the total coefficient of the term a.F.a 

3. K 10 = as it is the total coefficient of the term Tr[u.F.oj] 
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• Once K 8 , Kg and Ki are zero there are no terms in the fourth order divergence 
which are proportional to [u> 2 ] 2 , [a 2 ] 2 , a 2 u 2 or [a. a;] 2 . This will imply that in 
the divergence the net coefficients of all the terms, linear in a 2 , u 2 , (jj ilia uj av or 
a At u/ tly should be zero. 

To satisfy this condition we have to set all the from i = 1, • • • , 5 to zero. 
This will also set Qn to zero. 

Kq and Kf get related to B\ and B§ in the following way (see (B.4) and (B.5)). 

¥l = ¥l = l (^1 + ^1 + 2 dBl \ 
7] C s\dTTdTj 

Absence of these four fourth order terms mentioned above will also impose 
some constraints on the transport coefficients of second order stress-tensor by 
requiring that the coefficients of the four terms 0a 2 , 0a; 2 , oj a ba b uj ca and a a <J ab ab 
in the third order divergence should vanish. 



5 Constraints on 2nd order transport coefficients 

In this section we shall finally analyse how this condition of local entropy production 
constrains the second order transport coefficients. In the first part of this section we 
shall derive these constraints. These include the set of five relations among the 15 
transport coefficients (as mentioned in the section 1) and also two inequalities involv- 
ing both the first and second order transport coefficients as well as some coefficients 
appearing in the entropy current. 

Then in the next subsection we shall compare our final result with the answer 
presented in [2] and [1]. 



5.1 Derivation of the constraints 

At second order, just from symmetry analysis, the stress tensor will have 15 transport 
coefficients. 



n 



T 



t (u.V)a^ u) + KxR^v) + K 2 F{ p , v ) + A 0<x 



ij,i> 



+ Ai cr^ L a (x ai/ ) + A 2 a^ L a uj au ) + A 3 u^ a u au ) + A 4 a^a^ 



(5-1) 



+ TP. 



fJLV 



Ci(u.V)e + C 2J R + CsRoo + £i© 2 + 6<r 2 + ^u 2 + e 4 a 2 



As explained before, in the expression of the divergence of the entropy current, IP 1 ' 
will always appear contracted with a^ u and II with 0. Therefore, in U ab all the 
terms, which have either a a b or as factors, will finally generate a set of quadratic 
and higher order terms in a a b and 0. Such terms are always suppressed in dervative 
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expansion over the second order piece of the divergence provided the shear and the 
bulk viscosities are non zero. Therefore the coefficients multiplying these terms 
can never be constrained from the condition of positivity. Among the 15 transport 
coefficients, five (A , Ai, A 2 , Ci an d £2) are °f sucn type and therefore are completely 
unconstrained. 

It turns out that to maintain the positivity of the divergence, the coefficients 
r and £1 have to satisfy some inequalities. This is because at fourth order, the 
divergence of the entropy current will contain terms proportional to [(ti.V)a] 2 and 
[(m.V)©] 2 whose coefficients are Q 2 and Qi respectively (see (4.14)). These two 
terms, along with four other terms (a 2 , 2 , a tlu (u.'V)(T flu and 0(u.V)0, appearing in 
the second and third order pieces of the divergence) together can be made positive 
definite provided the transport coefficients r and Ci satisfy the following inequalities. 



(Ci - C e ) 2 < 4CQi 
(r - C CT ) 2 < 47yQ 2 



(5.2) 



Where Cq and C a are the coefficients of the term 0(w.V)0 and a ab {u.V)a ab respec- 
tively in the divergence of the third order entropy current. 



C = 2s— T-— + 2B 2 + 2£? 4 s s — T 



ds 3 dT \ dT 



(5.3) 



But unlike the inequalities for the first order transport coefficients (77 > and ( > 0) 
(5.2) involves several free coefficients appearing in the entropy current and hence it 
does not give any relation within the transport coefficients themselves. 

Now we shall come to those relations which will give some equalities among 
the remaining eight transport coefficients. By computing the divergence of the en- 
tropy current upto fourth order we can see that there are no terms proportional to 
R 2 , Rl , R ab R ab , F ab F a \ a 4 , u A and (a.tu) 2 . It will imply that the coefficients of the 
following 8 terms in the divergence of the entropy current have to be zero. 

1. C F = Coefficient of the term a ab F ab 

2. Cr = Coefficient of the term a ab R ab 

3. C a = Coefficient of the term a a a b a ab 

4. 0^ = Coefficient of the term u ap u p b a ab 

5. Qf = Coefficient of the term 0i?oo 

6. Qr = Coefficient of the term QR 
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7. Q a 

8. Q u 



Coefficient of the term Qa 2 
Coefficient of the term Quo 2 



Solving each of the above eight conditions we can express the remaining eight 
transport coefficients in terms of the coefficients appearing in the divergence of the 
second order entropy current. 

These are given by the following expressions. 



C R 

c,., 








A 3 = 



T 



A 3 , 
dB 5 



C F = 0^ k 2 = T 



dB 5 
dT 



dT 



4Bi 



C n = =>• Ad 



T 



,d 2 B, 



Qr 
Qf 

Qco 



o 
o 
o 



C 2 

£3 



dT 2 

dA 3 A 3 
ds 3 



T- 



^ + 2B 4 T 2 
dT 4 \dT 



dA 3 
ds 



3 

ds 



2T dB 5 



2B 4 Ts 



ds 
dT 



dB s 

-2B 4 Ts— + T— — ^ 
dT dT 



s dT 
T~ds 



dB l 
ds 



(5.4) 



Qa = => £4 



+ #1 



T 2 s 



2s dT 
T ds 

ds dB 4 



B A 



+ - [T 



where 



3 
dB 5 
dT 



dB, 



5 +T 2 



n 

3 

d 2 B 5 



ds \ 
dTJ 



4rr , ds ~ d 2 s 
+ 4Ts— + 2T 2 s— 
dT dT 2 



dT ' ' dT 2 
A 3 dA 3 
T dT 



From (5.4) we can see that all these eight coefficients can be determined in terms 
of three independent coefficients(A3, B\ and B 4 ) appearing in the third order entropy 
current. Therefore eleminating the three entropy current coefficients, mentioned 
above finally we shall get five relations among these eight transport coefficients. 
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K 2 

Cs 

6 



(is 3 



3 /s\ /dT 

4 vr) \~ds~ 



+ 



ck 2 



ds 
3 



3 

s 

f 



+ 2k 2 ) - ^ 
dT 2 4 



1 

+ 4 



A4 s / dT\ / T dA4 
"6 T \~ds ) [ X4+ 2dT 
Ts /dT\ fd 2 K 2 



T 



2 \ds J 



(d 2 K 2 \ 

\dT 2 J 



-) A 4 

ds J 



dn 2 

It 



Ts) Xa 



3s dT 
2T~ds~ 



(5.5) 



5.2 Comparison with [2] 

In [2] authors have constructed some examples of non-conformal fluid which can 
be obtained by dimensional reduction of some higher dimensional conformal theory. 
The entropy of such non conformal fluid is proportional to T 2u ~ x where 2a was the 
dimension of the space-time before the reduction. Since this particular nonconformal 
fluid satisfies the condition of 'positivity' of the divergence of the entropy current by 
construction the transport coefficients should also obey the relations listed in (5.5). 

Below we are quoting the values of some transport coefficients for such non- 
conformal fluids. These are the transport coefficients which enter the 5 relations in 
(5.5). 

A3 = A 3 T 2fJ 3 , £ 3 = 0/0 _ — — A 3 
Kl = kT 2ct " 3 , C2 = 77777 =rrKi 



a 4 = 0, e 



•1 



2(7- 


-4 


3(2a 


-1) 


2(7- 


-4 


3(2ct 


-1) 


2(7- 


-4 


3(2(7 


-2) 








(5.6) 

K 2 = (2(7 - 2)« 1 , C3 = ^777 777^2 



where A 3 and k are two dimensionful constants which depend on the length of the 
compactified dimensions but independent of temperature. Using the fact that for 
such dimensionally reduced nonconformal fluids the entropy can be written as 

rr2a— 1 
S OC T 

one can check that these values satisfy the relations given in (5.5). 
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5.3 Comparison with [1] 

To compare, first we shall express the eight relevant transport coefficients (the ones 
which appear in equation (5.5)) in terms of the coefficients as given in [1] . The 
dictionary is the following. 



TC2 


(-Rom r-p/- 

— ?5 ' 1 <■>: 




= -X Rom , 1 




T 2 (ds\ 


" s 2 [dTJ 



_ cRom rp(- tRom 
3 — ^6 ) 1 ^3 — — ^3 

n « 1 = K JtoB I TK 2 = 2(K-K*f 0m 

2 T 2 ( ds\ 2 

Rom rp\ _ ± I uo \ \Rom 



The author of [1] has argued for the existence of two relations among 5 of these 8 
nondissipative transport coefficients. These two relations are not explicitly presented 
in the paper [1] in generality. However the author of [1] appears to claim, in the un- 
numbered equation in section 5 (below equation 32 ) of [1], that in the special case, 
when 

T = s c % and K Rom ot^ (5.8) 
the two relations among the five transport coefficients reduce to the folllowing. 

^.Rom 

/3c 2 -1\ ^ 

cRom 1 cRom _ I s \ [ Rom . 2 \Rom~\ 

?6 + ?3 - _ I 3c 2 J L K + C s A 3 J 

The first equation in (5.9) indeed reduces to the second equation in (5.5) when c 2 
is a constant. In order to compare our results with the second of (5.9), we subtract 
the fourth equation from the third equation of (5.5) and then use the first equation 
of (5.5). 

This gives a relationship between all the same transport coefficients that appear 
in the second equation of (5.9). However the relationship we find is the following. 



Ca-6> 



ttiom , e. 
dK\ K\ 

3~ 



ds 



ds 



where 



k 2 = , 1 + T— 



1 / d«2 
3 VT 



3 s dT 
2TTs K2 



(5.10) 
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This relationship does not reduce to the second of (5.9) after substituing the special 
case of (5.8) with constant c s . We do not understand the reason for this disagreement. 
Perhaps the second of (5.9) applies under more restrictive assumptions than stated 
explicitly in [1]. As noted in [1] it certainly applies to the particular case, described 
in [2]. 

6 Conformal limit 

Upto second order in derivative expansion the final entropy current (consistent with 
the constraint of non-negative divergence) is given by the following expression 



J" 



second order 



=V„ [ii(«TT - u"VT)] + V, (A 2 Tu} ftu ) 

+ A 3 (V - l -g^R>j u u + + ^ [QWT - P Qb (V^)(V a T)] (6-1) 
+ (B x u 2 + B 2 Q 2 + B 3 a 2 )u" + B 4 [(Vs) V + 2s6V M s] 

If the theory has conformal symmetry, then the entropy current also should 
transform covariantly under a conformal transformation. The conformally covariant 
entropy current is a special case of equation (6.1). In this case the only available 
length scale is provided by the temperature and therefore the temperature depen- 
dence of all the coefficients are fixed just by dimensional argument and also some of 
the coefficients are related to the others in a way so that the terms that transform 
in-homogeneously under conformal transformation cancel. 

At second order in derivative expansion, there are three scalars and two vectors 
[1, 5, 8] which transform covariantly under conformal transformation. In our basis 
these are given by the following combinations 

= P ab V a V b T _ P ab (V a T)(V b T) _Rm_R & (6 2) 

3 j, 2T 2 2 4 6 

A conformally covariant entropy current should be expressible only in terms of 
these three scalars and two vectors. So to begin with it can have five independent 
coefficients. Then the constraint of positivity will reduce it to some special case of 
(6.1). 

Here we have used (6.1) to deduce the conformally covariant form of the entropy 
current. First we have fixed the temperature dependence of the coefficients Ai and 
Bi by dimensional analysis. Then we have tried to figure out the minimal set of 
relations these coefficients have to satisfy such that all the terms transforming in- 
homogeneously under conformal transformation cancel. This means that one should 



- 35 - 



be able to choose the coefficients Ai and Bi in such a way so that the entropy current 
is expressible in terms of these 3 conformal scalars and 2 conformal vectors. To do 
this we first rearrange some of the terms appearing in equation (6.1) assuming that 
the temperature dependence of the coefficients are fixed by dimensional analysis. 



A, 



^ - 2 (K + k, 



e 2 + #00 + 



R 



a b (a b » + cu b ») + |o" - \u k R ka P» a 



V, [A 2 uT\ = A 2 [V v 2 - S 2 u»] 

A 3 (et - \gTR) u v = A 3 -u» (| + #00) + P^R* 



(6.3) 



(6.4) 



T 



dAz 
dT 



[0V M T-P ab (V b ^)(V a T 



2A, 



e 5 



26 



-u 



(6.5) 



B A [(Vs) V + 2s9V M s] = T 6 B± 



9 V + 9 aV - 



20 



(6.6) 



From these expressions we can see how one should choose the coefficients A\ 
and Bi such that all the pieces that transform inhomogeneously under conformal 
transformation cancel. The coefficients for a conformally covariant entropy current 
are given by the following expressions. 



Ai(T) — ai, A 2 (T) = a 2 , A 3 (T) = ^-T 
£i(T) = bxT, B 2 (T) = ^T, B 3 (T) = b 3 T 



(6.7) 



B A {T) 



where all and 6, are constants. 

Therefore the conformally covariant entropy current has four independent coefficients 
(ai, a 2 , b\ and b 3 ) when expanded upto second order in derivatives. When written 
in terms of these four coefficients the expressions for the conformal entropy current 
is given as 



7^ 

conformal 



a{rs 3 u h 



aiT 



(Vf + V 2 M ) + a 2 T [V u 2 - S 2 u u ] + biTS 2 u" + b 2 TS lU ' 



(6.* 



T [ ai S 3 + b 2 S x + (h - a 2 ) S 2 ] u" + T (a 2 + ^) V£ + ^-V( 
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This expression coincides with the expression presented in [8] and [1] with the fol- 
lowing identification. 

a x T = 4A Rorn 



b 2 T 



3 

j^Rom J^Rom j^Rom 



2 

T{b x - a 2 ) = Af m + 2Af om - B* om (6.9) 



( j Rom 

2 

where Af om and Bf om are the coefficients in the conformal entropy current as defined 
in [1]. 

Substituting the relations (6.7) in (5.4) one can see that in conformal case £2, £3, 
£ 3 , £ 4 and A 4 vanish and k 2 is related to Ki as 

K 2 = 2k\ 

However once the stress tensor is conformally covariant, all these vanishing of 
the coefficients and the relation between k x and k 2 are automatic (If these relations 
were not true then the stress tensor would have some terms which will transform 
in-homogeneously under conformal transformation). Therefore we can say that the 
existence of an entropy with positive divergence does not constrain the uncharged 
conformal fluid. 
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8 Appendices 
A Identities 

Here we list the identities that we have used to calculate the divergence and then to 
transform the answer to the required basis. 

(u . v)s+se = 5^!±<5! + ... 

p,- v , r+ra ^v^+C9P-] + „, (AA) 
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26 



<x^(V M <v) = a*" [a°a au + u)^ a u av + + (u.V)<> - o M aJ + —a 2 (A.9) 



The identity (A. 10) is derived using the following steps. 

= uj^ u V u ^2c^ Q - u^a a + u a a^ 
= 2w^V i ,w Ma - tu 2 a a + (a^J u u a + w a V„a M ) 



[US 



'Rnaii.v R 



papv A ^ppau 



UJ 



V a u) vll + a^ u W a u u - — [R 



puap i - 1 ^ppa 



Rpaau] 



(A.10) 



V M CJ^A + ^^p\u + A Xufl 

- u uX a u - 2u u ^a x - (u p R pb P b u ) (a„6) 



where = u p R pabc P£P b u P x 



(A.ll) 



The identity (A.ll) can be derived using the similar tricks as in the identity 
(A.10). 

A^ uX V u [V x up + Vpux] 



A» uX V u 



2v\n + - u\<*p ~ u p a x 



A^ v ( - [V„, V A ] % + [V„ V„] u x + V M (V„« A ; 



(A.12) 



A pXu + A AjUI , + VpUv\ - QaV^m^ 



- 39 - 



B Computation of the divergence 



Here we shall calculate the divergence of the different terms appearing in the second 
order entropy current. The final expression for the second order entropy current is 
given in (6.1). Here we are quoting the equation again. 

J^ 1 1 second order 

=V„ [Ai(/VT - /VT)] + V u (A 2 Tu^ u ) 



±A 3 [R^--g^R)u v 



(B^ 2 + B 2 B< 



B 3 a 2 )u» 



y + iy) ^ v " t ~ pab (y^)(y- T f (R1) 

- B A [(Vs) V + 2s6Vs] 



The first two terms (with coefficients Ax and A 2 respectively) have zero divergence. 
Below we shall calculate the divergence of the rest of the terms. As explained be- 
fore, to determine the constraints of 'positivity' we need to calculate the divergence 
upto fourth order in derivative expansion. However in the fourth order piece of the 
divergence we need to retain only those terms which do not involve any factor of 
or 9. 

Divergence of the term with coefficient A 3 : 



upto 3rd order 



BT - ±<rR\ u u + A 3 (rT - \<TR 



{y,A 3 )R^u u - *(u.V)A 3 + A 3 (rT - 



6 



dA 3 A 3 \ R ( dA 3 

Roo{s ^7 + Y + 2 [ S Y7 



As 
3 



— l A* + T- 



dT 
V 



(uWR^) + A 3 R^ 



A 3 ( R>» - -gK"R ) u v 



- relevant terms at 4th order 



(B.2) 



Divergence of the term with coefficient (^f + 



v M [ev^ 5 -p a6 (v 6 ^)(v a B 5 

3 



ds 



(w.v)e 



a 



upto 3rd order 

dB 5 



T- 



dT 



a 2 6 



20 



e 



a^VaVvB* + — V 2 B 5 - —(u.V) Bi 



(B.3) 
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where 



dB 5 
dT 



A 3 dA 3 
T dT 



The three terms in the second line of equation (B.3) can be simplified further. 



(u.V) 2 B 5 



,d 2 B^ dB§ 



ds 2 



+ s- 



ds 



e 2 



dBx. . . „ „dBz. 



V 2 £? 5 



-e 2 



,d?B 5 TdB 5 



ds 2 



+ 



3 dT 



+ 



dB 5 
ds 



T- 



JB 5 
dT 



^B 5 
T—^ + T 
dT 



' a 2 - T^r(a 2 + oo 2 + iV) 



cff 2 



,dB* _ 2 c/ 2 5 5 



(IB* 2T dBs 



5 \ Ji 



ds 3 dT 



a 2 



2T— - + T 

The relevant part of the fourth order piece in the divergence is given by the 
following expression. 



v M [ewB 5 -p ab (v b u^(v a B t 



relevant part at 4th order 



dTh 
dT 



(B.4) 
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Divergence of the term with coefficient Bi. 



^ ' t 1 [^1^ ^ ] upto 3rd order 

B x u 2 Q + [{u^B^uj 2 - 2Bi 



, , J>, ,ca i ,,2® S dT 2 

U n h(J n U +u — — -co fc) 

3 2T ds 



dBi B 1 



ds 



a; 2 6 - AB^u^u, 



Ts 



2 /dBi 



s VdT 1 

7 \1t) + "t7 



1 1 v 



relevant part at 4th order 



r]P«V c a cb + CP ca V c 6 



25i 
7^ 



i^V^" - -i^V u e + R uu u u P^ 



pv" * a 



3' 

Divergence of the terms with coefficients B2 and S3: 

[5 2 eV] upto 3rd ordcr = e 3 (s 2 - ^+2^.v)e 



f 1 L ° J upto 



cr 2 e Lb 2 - s 



dB 2 



+ 2B 3 a flu (u.V)a 



I upto 3rd order ~ \ * ~ ds 

Divergence of the term with coefficient i? 4 : 

V M {B A [(Vs) 2 ^ + 2s6Vs]) I upto 3rd ordcr 

^ + B 4 ) 3 2 3 - 2fi 4 T 2 2 a,a u a^ + 2S 4 s6V 2 s 



- s 



+ a 2 6 



where V 2 s can be further simplified. 



V 2 s 



+ 



, ds\Q 2 ( ds\ . 

T ^lf)T + { a - T lf)^ 6 

d?s 1 a2 _ T & ^ + ^ 2 + ^ 



^ dS ^9 

T h T 2 - 



(B.5) 



(B.6) 



(B.7) 



There is no relevant part in the fourth order corrections to the equations (B.6) 
and (B.7) (i.e. all the terms appearing in the fourth order corrections to these 
equations involve atleast one factor of 9 or a^ u ). 
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